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The transport properties of interacting electrons for which the spin degree of freedom is taken 
into account are numerically studied for small two dimensional diffusive clusters. On-site electron- 
electron interactions tend to delocalize the electrons, while long-range interactions enhance localiza- 
tion. On careful examination of the transport properties, we reach the conclusion that it does not 
show a two dimensional metal insulator transition driven by interactions. A parallel magnetic field 
leads to enhanced resistivity, which saturates once the electrons become fully spin polarized. The 
strength of the magnetic field for which the resistivity saturates decreases as electron density goes 
down. Thus, the numerical calculations capture some of the features seen in recent experimental 
measurements of parallel magnetoconductance. 
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There has been much recent interest in the influence 
of electron-electron interaction (eei) on the localization 
properties of electrons in two dimensional disordered sys- 
tems. Behind this renewed interest in the topic are new 
experimental observations pertaining to the behavior of 
the conductance of low density two dimensional electrons. 
The conductance exhibits a crossover from an insulat- 
ing like temperature dependence. _at low densities to a 
metallic one at higher densities E'El A second transition 
back to an insulating dependence at even higher densi- 
ties was also observed.D This transition, which is known 
as the 2DMIT (two dimensional "metal-insulator" tran- 
sition), has drawn a flurry of theoretical activity since 
it is at odd with the prevailing single parameter scaling 
theory of localization.a This scaling theory asserts that 
for non-interacting electrons all states in 2D are local- 
ized by any amount of disorder. Since in low density 
systems the ratio between the typical interaction energy 
and the Fermi energy, rg, is large (i.e., > 1), a nat- 
ural explanation for the 2DM1T is that it is the result 
of the eei not taken into account in the original scaling 
theory. This has prompted an intensi\{e |theoretical ef- 
fort including analyticalcl and numerical3~Ej work which 
tried to explain the 2DMIT as a result of delocalization 
by the eei. On the other hand, one may argue that the 
observed temperature dependence of the conductance is 
not a result of a metallic zero temperature phase but 
rather a manifestation of essentially "high" temperature 
physics. Thus, some other physical mechanisiTa-|(|SYch as 
traps,til interband spin depepdeiit scattering,ta'E^ tem- 
perature dependent screeningElTtj or percolationE3) sets 
a very low temperature scale and the observed metal- 
lic behavior occurs at higher temperatures. Accordingly 
there is no 2DMIT and the systems are insulating at zero 



temperature. This viewpoint may find support in some 
recent experimental results which show a suppression of 
the 2DM1T once interband scattering is reduced,t3'E3 and 
from the observation that the bulk of the metallic behav- 
ior occurs at temperatures in which there is no quan- 
tum interference contributions to the conductance, while 
weak localizatioB^orrections are observed at very low 



temperatures .Eil~Eil 

Since it is extremely difficult to go beyond the per- 
turbative treatment of strong eei in disordered systems, 
many numerical studies have been performed in order to 
clarify the role played by the eei on the localization prop- 
erties of disordered systems. While for spinless electrons 
the eei qualitatively changes properties of r-tbe-, sjistem 
such as many particle energy^ eigl . st atistics ,E§Llj3'yj per- 
sistent current fiow pattern£3'LrEZIE3 pid charge density 
response to an external perturbation,t2l it does not lead 
to delocalization of electrons at the Fermi le\fel, j-which 
would manifest itself in enhanced conductivityBj^Eil 

The electron spin degree of freedom nevertheless plays 
an important role in the so called 2DM1T. When the 
influence of the eei on the conductance of disordered sys- 
tems is considered using a combinaticm. of perturbative 
and renormalization group techniques^ it leads to the 
conclusion that there is a divergence in the Cooperon 
channel at medium interactions, while the ladder chan- 
nel monotonously decreases as function of the interac- 
tion strength. Thus, for polarized (spinless) electrons 
eei should indeed decrease the conductance at zero tem- 
perature, while the situation for unpolarized electrons is 
unclear. 

The purpose of this paper is to investigate the influence 
of the eei on the transport through a disordered system 
once the spin degree of freedom is taken into account. In 
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previous studies in which spin was considered the effect 
of eei on various properties-of mesoscopic systems, such 
as the persistent currentcJu and the addition spectrum 
and spin polariza.tion of quantum dots in the Coulomb 
blockade regimeEjE^ was studied. In this paper we shall 
study the role played by eei on the tunneling amplitude. 
By evaluating the statistical properties of the tunneling 
amplitude of an interacting system one obtains informa- 
tion on the localization of the electrons, akin to the prop- 
erties revealed by the statistics of the wave functions of 
non-interacting elactrons, which is directly connected to 
the conductanceHj 

The effect of eei on the localization of the system is 
quite subtle. While the long-range part of the interac- 
tion always enhances the localization of the electrons, 
the short-range (Hubbard) part tends to delocalize the 
electrons. Thus, when both interactions are taken into 
account there is a competing influence, resulting in some 
delocalization for weak interactions while the localization 
is enhanced for stronger interactions. Although there are 
some superficial similarities to the 2DMIT, there are nev- 
ertheless some important differences, which in our opin- 
ion rule out delocalization due to the Hubbard interaction 
as an explanation for the 2DMIT. In order for a substan- 
tial delocalization to occur the Hubbard interaction must 
be significantly stronger than the long-range Coulomb 
part, i.e., the system is in the high density regime where 
good screening occurs. Indeed, for any reasonable esti- 
mation of interaction parameters, this delocalization oc- 
curs for electron densities corresponding to Ts < 2, while 
the 2DMIT occurs at much lower densities corresponding 
to rs ~ 4 - 20. 

Another striking feature exhibited by some systems 
for which a 2DMIT was observed is the strong decrease 
in the caruiiictance once an in-plain magnetic field is 
applied. eZlta This decrease coincides with the appear- 
ance of magnetizationEjO and saturates around the 
point for which the system becomes fully spin polar- 
ized. The saturation field decreases as the density goes 
down and becpmes zero at the "metal-insulator" tran- 
sition point .e1iE3 Naturally, one would like to relate the 
decrease in conductance to the appearance of spin po- 
larization. In fact, we find that the tunneling amplitude 
depends on the ground-state spin state once Hubbard in- 
teractions are considered, and since the application of an 
extecaal in-plain magnetic field changes the ground-state 
spirCJ it leads to a dependence of the conductance on the 
spin polarization. We observe a positive magnetoresis- 
tance, which is akin in some aspects to the experimental 
one, as well as saturation of the resistance at a critical 
magnetic field corresponding to full polarization. Unlike 
the enhancement in the conductance due to the Hub- 
bard interaction this behavior is robust to the addition 
of long-range Coulomb interactions. Moreover, the criti- 
cal magnetic field is reduced as the interaction strength 
is enhanced (i.e., corresponding to lower densities). Nev- 
ertheless, as we shall discuss later, a direct comparison 
between the numerical and experimental behavior leaves 



some open questions. 

We consider the following tight-binding Hamiltonian: 

[4j + l;^afcj;<T + 4+l,i;<Tafcj;^ + h.c] 

+U ^ {nk,j-a - K){ni^p^c,i - K)s/\rk,j - ri^p\, (1) 

where f = {k, j) denotes a lattice site, ^.^ is an electron 
creation operator (with spin a = —^,-1-^), the number 
operator is nkj-a — j cr^fej;^; ^k,j is the site energy, 
chosen randomly between —W/2 and W/2 with uniform 
probability, = 1 is a constant hopping matrix element, 
if = 1/ is a positive background charge equal to the elec- 
tronic filling and s is the lattice constant. As discussed 
above, there are two types of electron-electron interac- 
tion. The Hubbard interaction Uu due to the repulsion 
of electrons with opposite spin on the same site has a 
short range only, while the Coulomb interaction U has a 
long range. We include an in-plain magnetic field B, that 
couples only to the total z-component of the electron spin 

We consider systems composed of iV = 4 electrons re- 
siding on 6 X 6 lattices and systems of = 6 electrons 
on 4 X 4 lattices, corresponding to fillings oi i/ ^ 1/9 and 
3/8, respectively. We chose W = 5 and 8, respectively, so 
that the single particle localization length is comparable 
to the system sizes. Hard wall boundary conditions are 
chosen, since the lowest single-electron state for periodic 
boundary conditions is actually much more delocalized 
than neighboring states, resulting in an unusual behgu^ 
ior of the conductance close to the fully polarized stateEJ 
It is of course not possible here to directly mimic the ex- 
perimental procedure in which electron density is varied. 
Instead, the physical content of this density variation can 
be captured by controlling the ratio of the Fermi energy 
to the interaction energy. In the present model, it is 
achieved simply by changing the interaction strengths U 
and C/h while keeping other parameters constant. The 
value of U may be related to the electronic density via 
U = V^J Airvrg. On the other hand, there is no gener- 
ally accepted value for the ratio between U-r and U for 
2DEG. In Hubbard's original work the ratio was esti- 
mated as [/h = (10/3) {/-for weakly overlapping hydro- 
gen like wave-functions,L£l but there is no estimation of 
the ratio for relevant Si or GaAs samples. We there- 
fore will investigate both physical limits, U/Uh = 1 and 
U/Uyi = 0, as well as other intermediate values. 

We carry out our exact diagonalization in the subspace 
of the total number of electrons A^ and the total spin 
component Sz — M — N/2 with M being the number 
of spin up electrons. Since there is no mechanism for 



spin flip in the model, the many-particle wave functions 
with different values of Sz do not interact, and they can 
be calculated separately by block diagonalization. Us- 
ing the Lanczos method we obtain the many-particle 
eigenvalues e^'"^^ and eigenfunctions \a'^'^'). Because 
of spin symmetry, as long as -B = 0, e^'~^' = e^'^'' and 



for a given spin sector is (0^ 



The zero temperature local tunneling amplitude 
(0^|at^|0^^i) between the ground state of and A^- 1 
electrons can be employed here in order to characterize 
the transport properties of the many-particle interacting 
system. It has the advantage that only the ground state 
energy and eigenvector for N and — 1 electrons need 
to be calculated. The use of the tunneling amplitude in 
this context has-[ieen motivated and substantiated in our 
previous workMO If one compares the tunneling den- 
sity of states (TDOS) ^{e) in the independent particle 
approximation on the one hand and for the many-body 
interacting system on the other hand, it becomes evident 
that the tunneling amplitude (0^|at^|0^~^) replaces the 
single electron wave function. The same is true for the 
transmission t{f,r ',£,a) of an electron with energy e 
and spin a between two points r, r ' on the interface 
of the system with external leads, which is related to 
the conductance (T(e) through the Landauer formulatJ 
(j{e) = (e^//i) X^r.r ',o- i£i<^)\'^, wherethesumr, f ' 

is over all points on the interface. This behavior suggests 
that the tunneling amplitude is the appropriate quantity 
to replace the single electron wave function in studying 
transport properties of interaoting systems. A similar 
procedure is employed in Ref.Ej in order to generalize 
the concept of inverse participation ratio for interacting 
systems. 

Note, however, that once interactions are present a 
many particle state is a superposition of many differ- 
ent Slater determinants. Hence, the tunneling am- 
plitude is not normalized, as is the single electron 
wave function, which is the result of the fact that 
the spectral weight for interacting systems is not nec- 
essarily equal to one. Therefore, in order to study 
the influence of eei on quantum localization it is use- 
ful to define an effective tunneling amplitude 0(r, tr) = 
(O^lat jO^"i)/(Ef(0^l4,jO~-i)2)i/2. In the follow- 

ing, the effective tunneling amplitude (j){r, a) is traced 
for several Hubbard and Coulomb interaction strengths 
[/h and U as well as for different in-plain magnetic fields 
B. In order to analyze the tunneling amplitude and to 
derive the degree of localization, we calculate the inverse 
participation ratio P = „\4>{r)\'^ averaged over 100 
realizations. 

It is also possible to discuss the tunneling amplitudes 
between the lowest eigenvectors of a given spin sector 
Sz- Of course, when adding an additional electron, 
the total spin component Sz can either increase or de- 
crease by one half. Denoting the lowest eigenvalue for 
a given spin sector Sz with N electrons as |0^''^'), the 
tunneling amplitudes between the lowest eigenvectors 



with 5^"^' = Sf'^'""^ ± 1/2 for aU values of 4°"'^', i.e. 
^.initial = N ^ 4, &nd Sf'^'^""^ = for 

= 6. The results for the corresponding participation 
ratios are show in Fig. 1 as a function of C/h without 
Coulomb interaction [U = 0). While the tunneling am- 
plitudes of all possible spin channels have very similar 
participation ratios for C/h = 0, some of them are sig- 
nificantly increasing for larger Hubbard interaction, in- 
dicating weaker localization. In the case of completely 
polarized electrons, S'^ = | ^ 2 in (a) and S'z = | — > 3 
in (b), the Hubbard interaction has no effect, since it 
couples electrons of opposite spin only. It is interesting 
to note that there seems to be two distinct magnitudes 
of enhancement. The inverse participation ratio for tran- 
sitions in which S'f"'"-' = S'^'^"''^^ - 1/2 are substantially 
larger than for transitions where S'f"''' = S'™"''' -I- 1/2. 
This may be the result of the fact that for transitions with 
^.filial ^ 5'initiai -1 12 the additional spin-down electron 
can join many spin- up "partners", thus enhancing the 
effectiveness of the Hubbard interaction. For transitions 
with S*^"^! = S'i"'"''' + 1/2 the influence of the Hubbard 
interaction is diminished, because the relative number of 
possible pairs of electrons with opposite spin, that can 
share sites, is reduced. One may also observe that {P~^) 
for the 4x4 lattice shows some decrease for U-r > 5V, 
while for the 6x6 lattice more or less saturates. 

We attribute this behavior to the fact that for the smaller 
lattice the electronic density is higher (i^ — 3/8) resulting 
in an observable effect of the Mott-Hubbard transition, 
while for the lower density {v = 1/9) pertaining to the 
larger lattice the Mott-Hubbard transition is less pro- 
nounced. 

In order to find out which of the different possible 
spin channels is relevant for the zero temperature trans- 

port, the energies £a and £a ot the cor- 

responding multi-particle eigenstates have to be investi- 
gated. 5"^"'"'"' and S'f"''' have to be selected such, that 
both energies have their minimal value for each configu- 
ration. With no magnetic field, for a ground state spin 



value S, all energies 



-S-S+1,...,S-1,S 



are degen- 



erate. Once an infinitesimal magnetic field i? ^ is 
present (which we assume here), the ground state corre- 
sponds to the maximal value of 5*2 = S. Only if the 
condition = Si""''^^ ± 1/2 can be fulfilled after this 

minimizing procedure, the corresponding configuration 
has nonzero tunneling amplitude at zero temperature. 

Figure 2 shows the average values of the initial and 
final spins (open symbols, right scale) for zero tempera- 
ture transport as well as the average values of the partic- 
ipation ratio of the corresponding tunneling amplitudes 
(filled symbols, left scale) as function of the Hubbard 
interaction C/h, assuming B ^ 0^ . For the partici- 
pation ratios, three different averaging procedures are 
compared: arithmetic average {P~^), typical average 
exp — (InP), and geometrical average 1/{P). They give 
slightly different values, since the participation ratios are 



strongly fluctuating. But all of them for both considered 
systems show the same qualitative behavior: Upon in- 
creasing Hubbard interaction, the participation ratios for 
the zero-temperature transport are increasing as long as 
Uu < 5V, and they become weakly dependent on Uu for 
larger interactions. Hence, the zero-temperature trans- 
port in the 2DEG is enhanced by the Hubbard interac- 
tion. The enhancement is rather weak, though, reaching 
about 50% for the 6 x 6 lattice and 70% for the 4 x 4 
lattice. The range of Uu where the enhancement oc- 
curs, cannot be directly compared to the experiments, 
since the relevant parameter Vg is rather related to the 
Coulomb interaction U (which is zero here) than to the 
Hubbard interaction Uu. 

Next, we want to investigate, how the zero- 
temperature transport is affected by an in-plain magnetic 
field. While there would be some modification of the 
electeons' orbits for perpendicular magnetic fields (see 
e.g.Ell) the parallel magnetic field interacts only with the 
electrons' spins. Since the multi-particle eigenfunctions 
\a^'^') have already been calculated as eigenfunctions of 
the spin operator Sz, their energy levels e^''^^ are just 
shifted to 

- g^iBBSz according to Eq. (0). The 
eigenfunctions |a^''^^) and the corresponding tunneling 
amplitudes are not changes, but with B a differ- 
ent eigenstate might become the ground state. With in- 
creasing B, a higher degree of polarization corresponding 
to larger Sz becomes favorable.o Figures 3(a) and 3(b) 
show the average values of the ground state spins 5*™'*"^' 
(for the 3 electron system) and 5^ (for the 4 electron 
system) as function of the in-plain magnetic field B for 4 
values of Uu. In both figures, a monotonously increasing 
polarization is observed, as expected. 

Now, zero-temperature transport can be traced as 
function of B. One just has to average the participa- 
tion ratios of the tunneling amplitudes for the corre- 
sponding ground states, provided the condition Sf"*^' = 
S'™*'=^i±l/2 is fulfilled. If this condition fails, the tunnel- 
ing amplitude is zero for the corresponding configuration 
and B. The results of this i?-dependent averaging pro- 
cedure are shown in Fig. 3(c) for 3^4 electrons on 
a 6 X 6 lattice (upper part) and 5^6 electrons on a 
4x4 lattice (lower part), again for 4 values of Ub_. In the 
non-interacting case there is a weak dependence of the 
conductance on the magnetic field as a result of the weak 
dependence of {P~^) on the single particle state seen in 
Fig. 1 for f/n = 0. For nonzero Hubbard interaction Uu 
we observe a drop in (P~^) in the range of gfisB < A, 
where A is the mean single particle level spacing. This 
behavior is the result of the larger participation ratio ex- 
hibited by the S'f"^' = S'^'""''' - 1/2 transition compared 
to the Sf"''' = Sf"'*"' + 1/2 transition. As can be seen 
in Fig. 3, when g^BB < A/2, for all realizations, the 
ground state corresponds to 5*""'"^' = 1/2 for the odd 
ntimber of electrons, while more and more realizations 
change from S'f'^'^' = to S'f"''' = 1 as B increases. This 
leads to the decrease of the participation ratio of the av- 



erage tunneling amplitude, corresponding to a decrease 
in the conductance. The initial decrease is followed by 
intermediate maxima and minima, since, as B increases, 
the average Sz increases leading to a local peaks in the 
conductance each time that the ^'f'^^' = S'^"'"''' - 1/2 
becomes more prevalent. Once a larger portion of real- 
izations become fully polarized (for large magnetic fields 
B), the conductance drops to the same saturated value 
for all /7h values. The values Bsat of the magnetic field 
for which the systems become fully polarized are shown 
in Fig. 4(a) for both system sizes. It can be seen, that 
-Bsat decreases with increasing Hubbard interaction Uu, 
but it never vanishes. 

Our findings somewhat resemble the observations of re- 
cent experiments :Ej The conductivity of 2DEG semicon- 
ductor devices, that seem to show a "metallic phase", 
is reduced with increasing in-plain magnetic field, until 
a constant conductivity is reached for B > iJgat- Al- 
though this qualitative behavior is strikingly similar to 
the results of our calculation, the magnitude of the de- 
crease we observe is much smaller than that in the ex- 
periments. For silicon devices decreases by up to two or- 
ders of magni±j|id£ have been observed depending on the 
temperature,L3o while for GaAs,^vices, decreases by a 
factor of 3 have been reported.oEy Another problem in 
relating our results to the experiments is the dependence 
of the saturation field Bgat on Uu- In the experiments 
^sat c>c — , where z/c is the densitw,of the 2DMIT 
and (5 = 1 for a wida-^nge of densitiesj^d while S « 0.6 
close to Vc elsewhere.Ea Since we do not identify a metal- 
insulator transition in our model it is not clear how to 
directly relate the above experimental observations to our 
data. Moreover, the experimentally relevant parameter 
is rather related to the Coulomb interaction U (via r^) 
(which is zero here) than to the Hubbard interaction Uu- 

Nevertheless, it is possible to relate i?sat to Uh using 
the following consideration: The lowest average energy 
of the Hamiltonian giwen in Eq. |l| at a given spin S may 
be approximated by:c3 

E{S)~E{0)^AS^~JS{S + l)-gnBBS, (2) 

where J is the averaged exchange energy which depends 
on U and Ub- Figure 4(b) shows the values of J de- 
termined from the eigenvalues E{Sniax) and E{0) of the 
systems via Eq. (||) for i? = and U = 0. A power-law re- 
lation J (X U^ with an exponent a in the range from 0.25 
to 0.5 may be deduced from our data, and this relation 
is consistent with results obtained, in a previous study 
of the disordered Hubbard model.L3 The saturation field 
Ssat corresponds to the field in which the energy becomes 
minimal for complete polarization, i.e. dE{S)/dS = 
with S = S'rnax, resulting in 

g/^B^sat = 2A5max " J (2S',„ax + 1)- (3) 

For small values of S it is better to replace the deriva- 
tive leading to Eq. (||) by the discrete difference leading 
to 5/ZBSsat = A(25'max " 1) " aJS^ax- In Figure 4(b) 



the value of the averaged exchange energy J obtained 
from i?sat via this equation is compared to the J ob- 
tained from the eigenvalues. The Figure shows that the 
relations work reasonably well in describing i?sat- Since J 
increases as function of C/h the saturation field decreases. 
Thus, the same physics of complete spin polarizaticm-at 
the saturation field describes both the experimentEilO 
and the numerical results. It is important to note though, 
that since the exchangalj saturates at values lower than 
A, J{Uy{ oo) ~ 0.6A, there will be no ferromagnetic 
transition for any value of C/h- This is clear from Fig. 2 
where the average value of the spin saturates at low val- 
ues. Therefore i?sat 7^ for any value of C/h, and a model 
with no long-range interactions can not describe the ex- 
perimentally observed -Bgat oc {ly — I'cY . 

Thus, it is clear that one must move beyond the 
short-range Hubbard interaction and include a long- 
range Coulomb interaction in order to present a realistic 
picture of the 2DMIT behavior. As discussed above there 
is no established relation between U and C/h, so we shall 
begin by presenting the highest possible ratio between 
them, i.e., U = C/h- The participation ratio as function 
of U is presented in Fig. 5(a). An initial enhancement 
of P^^, which peaks at C/ = 2V is seen, while for larger 
values of interaction a steady decrease is observed. This 
behavior is attributed to the interplay between the short 
and long-range components of the eei. Short-range inter- 
actions enhance the participation ratio of the tunneling 
amplitudes, while the long-range interactions suppress 
it. This is demonstrated in Fig. 5(b) where we keep 
C/h — lOV fixed and tune the value of U. A monotonous 
decease in as function of U is evident. The influence 
of the long-range interaction on the spin polarization can 
also be gleaned from Fig. 5. The short-range interactions 
lead to spin polarization, while the long rang interaction 
tends only to slightly nudge the ground state spin. 

The dependence of the participation ratio on the in- 
plain magnetic field in the presence of long range eei is 
depicted in Fig. 6. The main feature of Fig. 3, i. e., 
the suppression of by the magnetic field and satu- 
ration at high magnetic fields remains intact. Thus, long 
range Coulomb interactions do not overturn the main 
conclusions garnered from the magnetoconductance for 
short-range interactions. The dependence of the satu- 
ration field Bsat on the Coulomb interaction strength U 
shown in Fig. 4(a) is also similar to the dependence on 
C/h for U — 0. There are some differences though. The 
saturated value of decreases as U becomes stronger, 
which fits our previous observation that long-range inter- 
action localizes the system. The saturation field depends 
on the interaction strength, but in contrast to the Hub- 
bard interaction we can not rule out that Sgat = for a 
finite value of C/. If Ssat = indeed occurred, the value 
of interaction for which the field is equal to zero would 
be much larger than the values of interaction considered 
here. Therefore, we cannot compare the behavior of J5sat 
to the one seen in the experiment based on the available 
data. 



Once long-range interactions are taken into account 
one can use the relation U = V^/AttiTs, resulting in a 
connection between interaction strength and electronic 
density. It would be tempting to interpret Fig. 5(a) as a 
signature of a 2DMIT, but this explanation raises diffi- 
culties which make it rather doubtful. The first "transi- 
tion" from a more insulating behavior to a more metallic 
one (see Fig. 5(a)) is influenced almost entirely by the 
Hubbard term as can be seen from its similarity to Fig. 
2. At C/ « 2V, the point where the enhancement due 
to the short-range interaction saturates, the suppression 
due to long-range interactions kicks in. Since the satura- 
tion point is determined by C/h, while depends on U 
the corresponding electronic density for which the peak 
appears to depend on the ratio of C/h/ C/. For the lowest 
possible ratio (i.e., Uu/U = 1) the peak occurs at densi- 
ties corresponding to ^ 1.6 which is low compared ta 
the experimentally observed region of metallic behaviond 
corresponding to rs ~ 4 — 20. Similar behavior for the 
site occupation number, which gives some indication on 
the degree of localip.tion in the system, has been seen by 
Selva and Pichard.Efl A higher ratio of C/h/C/ will result 
in this peak appearing for even lower values of r^. The 
high densities (low values of r^) for which the peak in 
the participation ratio appears, in conjunction with the 
fact that long-range Coulomb interactions always sup- 
press P~^ (see Fig. 5(b)) lead us to the conclusion that 
numerical studies of small clusters do not support the 
notion of an eei driven metal insulator transition. 

On the other hand, the numerical model does seem 
to reproduce the positive in-plain magnetoresistance and 
the saturation of the resistancepat_a critical magnetic 
field seen in recent experimentsJllTEj The physical ori- 
gin of the saturation is spin polarization of the elec- 
tron. Similar explanations were prpposcd for the ex- 
perimental origin of the saturation.EirEII Further stud- 
ies of higher values of interactions are needed to deter- 
mine whether quantitative comparison with the exper- 
imentally observed i?sat oc {ly — VcY can be obtained. 
Furthermore, if such behavior was indeed observed, the 
origin of Vc would have to be clarified, since it can not 
be the critical density of the 2DMIT which is not seen in 
our model. 

Beyond relating our results to existing experimental 
data, one would like to glean some relevance for possible 
future experiments. As we pointed out in the introduc- 
tion the physical situation corresponding to a short range 
eei is deep in the metallic regime in which the long-range 
interaction is perfectly screened. Thus, the enhancement 
of the conductance in Fig. 2 is not relevant to the ex- 
isting body of experimental data on the 2DMIT, but it 
might be relevant to an experimental double gate set-up 
in which the long-range part of the eei is screened.E3 Ac- 
cording to Fig. 2 the conductance of such a set-up should 
be higher than for a single or no gate set-up for the same 
2DEG sample. 

In conclusion, the conductance of interacting electrons 
for which the spin degree of freedom is taken into account 



exhibits an intricate dependence upon eei and parallel 
magnetic field. Hubbard on-site interactions enhance the 
conductance while long-range Coulomb interactions sup- 
press it. The interplay between the two can lead to a 
region of enhanced conductance, but this region is at den- 
sities corresponding to « 1, which are too high to be 
relevant to the 2DMIT. A parallel magnetic field reduces 
the conductance until it saturates once the electrons be- 
come fully spin polarized. The saturation field decreases 
for lower electron density. These features are seen in 
recent experimental measurements of parallel magneto- 
conductance. Further investigation, at higher values of 
interaction, are needed in order to clarify whether other 
features may be explained by this model. 
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UJV 

FIG. 1. The participation ratio P ^ of the tunneling amplitude for the different spin channels is shown as function of the 
Hubbard interaction strength Uh for (a) 3 ^ 4 electrons on a 6 x 6 lattice, and (b) 5^6 electrons on a 4 x 4 lattice, both 
without Coulomb interaction. The symbols correspond to the spin transitions: 5^ = | — > (squares), f — > 1 (circles), § ^ 1 
(triangles up), | ^ 2 (triangles down), | — > 2 (diamonds), | — > 3 (plus). We have averaged 100 realizations of disorder, and 
the error bars show the standard deviations of the averages. 
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FIG. 2 

average initial and final spin (open symbols, right axes) are shown as function of the Hubbard interaction strength Uu for 
(a) 3 — > 4 electrons on a 6 x 6 lattice, and (b) 5^6 electrons on a 4 x 4 lattice, both without Coulomb interaction. In the 
averaging procedure, the ground states for both, initial and final state have been chosen as described in the text. For all 100 
disorder configurations, the tunneling amplitude is non-zero. For the participation ratios, arithmetical, typical, and geometrical 
averages are compared. 




FIG. 3. The average (a) initial and (b) final spins Sz of the ground states are shown versus the scaled magnetic field for 3^4 
electrons on a 6 x 6 lattice. With increasing magnetic field, the ground states are moved to higher degrees of polarization. The 
different lines correspond to different strengths of Hubbard interaction Un (see legends), while there is no Coulomb interaction. 
On the right, the average paxticipation ratios (-P~^) of the zero temperature tunneling amplitude are shown for 3 — > 4 electrons 
on a 6 X 6 lattice (c) and 5 — > 6 electrons on a 4 x 4 lattice (d). In the (arithmetic) averaging procedure, the ground states for 
both, initial and final state have been chosen as described in the text. Configurations with zero tunneling amplitude (due to 
non-matching spins in initial and final ground states) have been disregarded in the averaging procedures. 
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FIC. 4. (a) The saturation field Bsat, for which the systems become fully polarized and the average participation ratios 
of the zero temperature tunneling amplitude reach their asymptotic value, are shown versus the Hubbard and the Coulomb 
interaction strength (see legend). The values of Bs&t have been averaged for 100 configurations, (b) The values of the average 
exchange energy J are shown versus ?7h for 6x6 systems (open symbols) as well as the 4x4 systems (filled symbols). Two 
ways to calculate J are compared: (i) J determined from the eigenvalues via Eq. (2) for B = (squares) and (ii) J obtained 
from Bsax via the discrete version of Eq. (3) (circles). The plot indicates good agreement, justifying the approach to relate Bgat 
to U-R. 
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FIG. 5. The average participation ratio P ^ of the zero temperature tunneling amplitude (filled symbols, left axes) and the 
average initial and final spin Sz (open symbols, right axes) are shown as function of Coulomb interaction strength U, (a) for 
identical Hubbard interaction Uh = U and (b) for fixed Hubbard interaction Un = WV for 3 ^ 4 electrons on a 6 x 6 lattice. 
The averaging procedure is the same as for Fig. 2. 
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FIG. 6. The average participation ratio {P ^} of the zero temperature tunneling amplitude is shown versus the magnetic 
field B, (a) for identical Hubbard interaction Uu — U and (b) for fixed Hubbard interaction Uu = 10 for 3^4 electrons on a 
6x6 lattice. The averaging procedure is the same as for Fig. 3. 



